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^ Abstract 



The interacting polytropic gas dark energy model is investigated from the view- 
ed 

point of statefinder diagnostic tool and w — w analysis. The dependency of the 
statefinder parameters on the parameter of the model as well as the interaction pa- 
rameter between dark matter and dark energy is calculated. We show that different 
values of the parameters of model and different values of interaction parameter result 
different evolutionary trajectories in s — r and w — w' planes. The polytropic gas 
model of dark energy mimics the standard ACDM model at the early time. 



* Email: malekjani@basu.ac.ir 
t Email: khodam@basu.ac.ir 



2 



I. INTRODUCTION 



Recent astronomical data obtained by SNe la [l], WMAP 2\, SDSS |3] and X-ray 4] 
experiments suggest that our universe expands under an accelerated expansion. In the 
framework of standard cosmology, a dark energy component with negative pressure is re- 
sponsible for this acceleration. A major puzzle of cosmology is the nature of dark energy 
and therefore many theoretical models have been proposed to interpret the behavior of dark 
energy. Although the earliest and simplest model is the Einstein's cosmological constant, 
but it suffers from tow deep theoretical problems namely the "fine-tuning" and "cosmic co- 
incidence" . The other models for dark energy scenario are the dynamical models in which 
the EoS parameter is time-varying. According to some analysis on the SNe la observational 



data, it has been shown that t 
with a cosmological constant 



re time- varying dark energy models give a better fit compare 
ij]. The dynamical dark energy models are classified in tow 



js], K-essence {9], tachyon [10], dilaton Jill 



energy such as Chaplygin gas models 12 



different categories: (i) The scalar fields including quintessence [6], phantom |7j, quintom 

and so forth, (m) The interacting models of dark 
13], braneworld models [yj], holographic [l5] and 







agegraphic 16] models. The interacting dark energy models have also been investigated in 
17]- 



The holographic dark energy model comes from the holographic principle of quantum 
gravity [18] and the agegraphic model has been proposed based on the uncertainty relation 
of quantum mechanics together with general relativity 19]. Recent observational data gath- 
ered from the Abell Cluster A586 support the interaction between dark matter and dark 
energy [20 1. However the strength of this interaction is not clearly identified 21]. 
Here in this work we focus on the polytropic gas model as a dark energy model to explain 
the cosmic acceleration. In stellar astrophysics, the polytropic gas model can explain the 
equation of state of degenerate white dwarfs, neutron stars and also the equation of state of 
main sequence stars 22j . The idea of dark energy with polytropic gas equation of state has 
been investigated by U. Mukhopadhyay and S. Ray in cosmology 23]. Recently, Karami 
et al. investigated the interaction between dark energy and dark matter in polytropic gas 
scenario, the phantom behavior of polytropic gas, reconstruction of f{T)- gravity from the 
polytropic gas and the correspondence between polytropic gas and agegraphic dark energy 
model [24j-l26| . The cosmological implications of polytropic gas dark energy model is also 



discussed in {27)]. The evolution of deceleration parameter in the context of polytropic gas 
dark energy model represents the decelerated expansion at the early universe and acceler- 
ated phase later. Depending on the parameters of the model, the polytropic gas can achieve 
a quintessence regime. The potential and the dynamics of K-essence, dilaton and tachyon 
fields according to the evolution of polytropic gas model is also investigated in {27 1. 
The polytropic gas is a phenomenological model of dark energy. In a phenomenological 
model, the pressure p is a function of energy density p, i.e., p = —p — f(p) j^. For 
f(p) = 0, the equation of state of phenomenological models can cross w = —1, i.e., the 
cosmological constant model. Nojiri, et al. investigated four types singularities for some 
illustrative examples of phenomenological models [28] . The polytropic gas model has a type 
III. singularity in which the singularity takes place at a characteristic scale factor a s . 
Since many theoretical dark energy models have been proposed to explain the accelerated 
expansion of the universe, therefore a sensitive test which can discriminate between these 
models is required. The Hubble parameter, H = a/ a, (first time derivative) and the de- 
celeration parameter q = —d/aH 2 (second time derivative) are the geometrical parameters 
to describe the expansion history of the universe. Since a > 0, hence H > means the 
expansion of the universe. Also a > 0, i.e. q < 0, indicates the accelerated expansion of 
the universe. Since the various dark energy models give H > 0, q < at the percent time, 
hence the Hubble parameter and deceleration parameter can not discriminate dark energy 
models. For this aim we need a higher order of time derivative of scale factor. Sahni et al. 



29[ and Alam et al. [30], by using the third time derivative of scale factor, introduced the 
statefinder pair {s,r} in order to remove the degeneracy of H and q at the present time. The 
statefinder pair has been given by 

a r-l 

? s = o^ I/on (!) 



aH 3 ' 3(g-l/2) 

Depending the statefinder diagnostic tool on the scale factor indicates that the statefinder 
parameters are geometrical. Up to now, the various dark energy models have been studied 
from the viewpoint of statefinder diagnostic. The various dark energy models have differ- 
ent evolutionary trajectories in {s, r} plane, therefore the statefinder tool can discriminate 
these models. The well known ACDM model is related to the fixed point {s=0,r=l} in the 
s — r plane [29I ]. The other dynamical dark energy models that have been investigated by 
statefinder diagnostic tool are: 
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the quintessence DE model 
graphic dark energy models 



29 
33 



3C| , the interacting quintessence models [3l|, |32|, the holo- 
34( , the holographic dark energy model in non-flat universe 



351 ] . the phantom model [36[, the tachyon [37]], the generalized chaplygin ga s model [38[, the 



interacting new agegraphic DE model in flat and non-flat universe 39 



401 ] . the agegraphic 



dark energy model with and without interaction in flat and non-flat universe 44 



421 ] are 



analyzed through the statefinder diagnostic tool. Recently, the statefinder parameters have 



been investigated by considering the variable gravitational constant G 



43]. 



In addition to statefinder diagnostic, the other analysis to discriminate between dark energy 



models is w — w' analysis that have been used widely in the papers 4CH42 



44 



5311 . In this 



work we investigate the interacting polytropic gas model by statefinder diagnostic tool and 
w — w' analysis. We introduce the interacting polytropic gas dark energy model in sect. II. 
The numerical results is presented in sect. III. We conclude in sect. IV. 

II. INTERACTING POLYTROPIC GAS DARK ENERGY MODEL 



In this section we give a brief review of the polytropic gas model for dark energy scenario. 
For more details and discussion see j^]. The equation of state (EoS) of polytropic gas is 
given by 



(2) 



22] 



where K and n are the polytropic constant and polytropic index, respectively 
In the non-flat FRW universe including dark energy and dark matter components, the Fried- 
mann equation is given by 



H 2 + - 2 

a z 



3M p 2 



(Pm + Pa) 



(3) 



where H is the Hubble parameter, M p is the reduced Planck mass and k = 1,0,-1 is a 
curvature parameter corresponding to the closed, flat and open universe, respectively. p m 
and pa are the energy densities of dark matter and dark energy, respectively. 
The dimensionless energy densities are defined as 

k 



Pc 



Pr, 



Pa 



3M p 2 # 2 ' " ,l p c 3M2H 2 
Therefore the Friedmann equation ([3]) can be written as 



a 2 H 2 



(4) 



(5) 
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In a universe dominated by interacting dark energy and dark matter, the total energy density, 
P = Pm + Pa, satisfies a conservation equation 

p + 3H(p + p) =0 (6) 

However, by considering the interaction between dark energy and dark matter, the energy 
density of dark energy and dark matter does not conserve separately and in this case the 
conservation equations for each component are given by 



p m + 3Hp m = Q, 
Pa + 3H(p k + p A ) = -Q, 



(7) 
(8) 



where Q represents the interaction between dark components and can be given as one of the 
following forms 54 ] 

3aHp A 

Q = I 3f3Hp m (9) 
^ 3'jH(p A + p m ) 

where a, /3 and 7 are the dimensionless constants. The Hubble parameter H in the Q-terms 
is considered for mathematical simplicity. The interaction parameter Q indicates the decay 
rate of dark energy to the dark matter component. The interaction between dark energy 



and dark matter is also studied in [55j. Here same as our previous worksjfor mathematical 
simplicity, we consider the first form of interaction parameter Q {27! . 40, 42]. Using Eq.(J2]), 
the integration of continuity equation for interacting dark energy component, i.e. Eq.(JH]), 
gives 



Pa 



1 



3(1+00 ~ 

Ba n — K 



(10) 



where B is the integration constant, K = and a is the scale factor. A positive energy 
density for any value of n is achieved when Ba 3 ^ 1+a ^ n > K. In the case of even n, we 
have positive energy density for any condition of K. The phantom behavior of interacting 
polytropic gas dark energy has been studied in 25J. In the case of Ba 3(1+a)/n = K, we have 
p — > 00 and the polytropic gas has a finite-time singularity at a c = (K / B) n ^ 3( - 1+a \ This 
type of singularity, in which at a characteristic scale factor a s , the energy density p — > 00 
and the pressure density \p\ — > 00, has been indicated by type III singularity 28]. 
Substituting Q = 3aHp A in (|Sj) obtains 



Pa + 3H(1 + a + w a )pa = 0, 



(11) 
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Taking the derivative of Eq. lflO]) with respect to time, one can obtain 

p A = -WH(l + a)a^^ ii pj v + " (12) 

Substituting Eq. ffl2|) in (TTTT) and using Eq. lflOl) , we can obtain the EoS parameter of inter- 
acting polytropic gas as 

3(l + a) 
Ob n 

w\ = -l -a (13) 

c — a n 

where c = K/B. We see that the interacting polytropic gas model behaves as a phantom 
model, i.e. w\ < —1, when c > a 3 ( 1+Q 0/ n . It is also clear to see that at the early time 
(a — y 0) and the absence of interaction between dark matter and dark energy (a = 0), the 
polytropic gas mimics the constant, i.e. w\ — > — 1. The evolution of EoS parameter can be 
obtained by differentiating of f fT3"j) as follows 

, 3(1 + a)ca n 

< = — Wk7 ( 14 ) 

n{c — a n y 

where prime denotes the derivative with respect to x — In a. Using Eqs.f llOp and @ the 
density parameter of interacting polytropic gas is given by 

= mpp (15) 

Taking the time derivative of Eq.f JT5|) and using Q' = Q/H, yields 

»^-»a( 3(1 ^^ + ^) (16) 
v a n — c n 7 

Taking the time derivative of Friedmann equation fl3]) and using Eqs. ffTOj) . ([5]), ((7|), f fT5|) and 

Q = 3aHp\, one can find that 



H_ 3 



c(l + a) n fc i 

"A 3(l+a) h "3~ 



(17) 



a « — c 

Substituting this relation into Eq.( Tl6|) . we obtain the evolutionary equation for energy den- 
sity parameter of interacting polytropic gas as: 

fi A = -3fi A - ^^y (1 - Ov) + « 3 (i +a) T" > ( 18 ) 

L a « — c a « — c J 



The evolution of density parameter has been discussed in [27|. It has been shown that 
the polytropic gas dark energy model can describe the matter dominated universe at the 



7 



early time, fi A — > 0, and dark energy universe at the late time, fi A — > 1, see Fig. (2) of [27] . 
The deceleration parameter is given by 

a H . . 



Substituting ( flTI) in ( fl9i) . the deceleration parameter can be obtained as 

(20) 



, 3 



c(l + a) ^ 

"A 3(l+q) ^ 1 + 

a n — c ° 



The evolution of q for interacting polytropic gas model has also been presented in [27]. At the 
early time the universe has a decelerated expansion, q > 0, and enters into the accelerated 
phase later, q < 0, see Fig. (3) of j^ . 
Using (H3D, ([IZD and ([20D, we have 

H 9 

— = --^a(1 + a) (a + wa)[(1 + a)(-w A + fi A a + fi A w A ) - a(a + 2)] 

-|fi A (l + (*Wa + |["a(1 +«)(« + tw A ) + I] 2 (21) 

At what follows, we derive the statefmder parameters (s,r) for polytropic gas model in 
the interacting spatially flat universe. Using the definition of statefmder parameters in (JT]), 
one can obtain 

U 3q-2 (22) 



aH 3 H 3 

Inserting fTSOj) and (j2ip in fT22|) and using fTl6|) we have 

3 

r = 1 + -n A (l + a) [3(1 + a) (a + to A )(l + a + w A ) - u; A ] (23) 

Inserting ( 1201) and ( |23l) in ([T]), the parameter s for interacting polytropic gas is obtained as 

_ 2 3a(a + l) 2 + 3aw A (2a + w A + 3) + 3iy A (l + w A ) - w' A 
3 a + w A 

In the limiting case of w A = —1, it is obvious w' A = and in the absence of interaction 

between dark matter and dark energy, i.e. a = 0, the statefmder parameters reduce to 

{s = 0, r = 1} which is coincide to the location of standard ACDM model in s — r plane. 



III. NUMERICAL RESULTS 



In the present section we give the numerically description of the evolutionary trajectories 
of the statefmder parameters in s — r plane for interacting polytropic gas dark energy 
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model in the flat universe. We also perform the w — w' analysis for this model. Here we set 
fl^ = 0.3 and Q° =0.7 for the density parameters of dark matter and dark energy at the 
present time, respectively. 

A. Statefinder diagnostic 

The statefinder pair {s, r} in this model is given by (12 3 p and §M§- One can easily see the 
dependency of the {s,r} on the EoS parameter , w, as well as the interaction parameter a 
in (123]) and (l2~4"j) . In the limiting case of non-interacting polytropic gas (a = 0), from ffTB"]) 
we see that at the early time (a — > 0) the EoS parameter % — > — 1. From ffl^j) . we also 
have w' A — > 0. From f )23|) and (j2l|) we see that at the early time the statefinder parameters 
for non-interacting flat universe are (s = 0,r = 1) which is coincide to the location of 
spatially flat ACDM model in the s — r plane. Hence, the polytropic gas model mimics the 
ACDM model at the early time. 

In Fig.(l), the evolutionary trajectories of interacting polytropic gas model is plotted 
for different values of interaction parameter a. Here we fix the parameters of the model 
as c = 2,n = 4. The standard ACDM fixed point is indicated by star symbol in this 
diagram. The colored circles on the curves show the present values of statefindr pair 
{sq,t }. Different values of a result different evolutionary trajectories in s — r plane. Hence 
the interaction parameter can influence on the evolutionary trajectory of polytropic gas 
model in s — r plane. For larger value of a, the present value So decreases and the present 
value ro increases. The distance of the point (so,To) form the ACDM fixed point (i.e. 
s — 0,r = 1) becomes larger by increasing the interaction parameter a. While the universe 
expands, the evolutionary trajectory of interacting polytropic gas dark energy model 
evolves from the ACDM at the early time, then r increases and s decreases. The present 
value {so, r } are valuable, if it can be extracted from the future data of SNAP (SuperNova 
Acceleration Probe) experiments. Therefore, the statefinder diagnostic tool with future 
SNAP observation are useful to discriminate between various dark energy models. 
In Fig. (2), the evolutionary trajectories for interacting polytropic gas are plotted for 
different values of the parameters of the model. Here we fix the interaction parameter as 
a = 0. In left panel, the parameter n is fixed and the parameter c is varied. Different values 
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of c gives the different evolutionary trajectories in s — r plane. Therefore the parameter c 
of the model can affect on the evolutionary trajectories in s — r plane. Like Fig.(l), the 
present value of statefmder pair, i.e. {so,r } is indicated by colore circles on the curves. 
For larger values of c, r$ decreases and so increases. The distance of the point (sq, tq) to the 
location of standard ACDM fixed point becomes shorter for larger value of c. In right panel 
the parameter c is fixed and the parameter n is varied. Same as left panel, the interaction 
parameter is fixed to a = 0. Here we also see that different values of n gives different 
evolutionary trajectories in s — r plane. For larger values of n, we see r decreases and So 
increases. Here we see that, same as parameter c, the distance of the point (so,ro) to the 
location of standard ACDM fixed point becomes shorter for larger value of n. 



B. w/\_ — w' A analysis 

In addition to the statefmder diagnostic, another analysis to discriminate various models 
of dark energy is w — w' A analysis. Here we apply this analysis for interacting polytropic 
gas dark energy model. In this analysis the standard ACDM model corresponds to {w\ = 
—l,w' A = 0}. The evolution of wa and w' A are given by (fT3l and (jLT| . respectively. In 
Fig. (3), the evolutionary trajectories of interacting polytropic gas dark energy for different 
values of interaction parameter are shown in w\ — w' A plane. Here we fix the parameter 
of the polytropic model as c = 2 and n = 4. One can see the different values of a result 
different trajectories in wa — w' A plane. The present value w° A — w A is dependent on the 
interaction parameter. Larger value of a obtains smaller values of w A and w A . 
In Fig. (4), the evolutionary trajectories are plotted in the absence of interaction parameter, 
i.e. a = 0. Here we perform the — w' A analysis for different values of the parameters of 
model. In left panel the parameter n is fixed and in the right panel the parameter c is fixed. 
In these diagrams, by expanding the universe, the evolutionary trajectories start from the 
fixed point wa = —l,w' A = (i.e. the location of ACDM fixed point). In left panel, one 
can see that the parameter c can only affect the present value w A — w A . Different values 
of c result the same evolutionary trajectory in w a — w' A plane. The distance of the present 
value w A — w' A to the location of ACDM fixed point (i.e. {wa = —^,w A = 0}) is shorter 
for larger value of c. In right panel we have different evolutionary trajectories in w\ — w' A 
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plane for different values of n. The parameter w'® increases for larger value of n. Like left 
panel, the distance of the present value w° — to the location of ACDM model is shorter 
for larger value of n. 



IV. CONCLUSION 



Since many dynamical dark energy models have been proposed to interpret the cosmic 
acceleration, the statefinder diagnostic tool based on third time derivative of scale factor 
is given to discriminate between them. The statefinder diagnostic combined with future 
SNAP observation can be useful to discriminate between various dark energy models. Here 
we studied the statefinder diagnostic tool for interacting polytropic gas model in spatially 
flat universe. We derive the statefinder parameters s and r in this model and investigate 
the dependency of the evolutionary trajectories in s — r plane on the parameters of the 
model as well as the interaction parameter between dark matter and dark energy. We 
obtained the present value {so, r } of this model and studied the dependency of {so, r } 
on the parameters of the model and interaction parameter. By expanding the universe, the 
evolutionary trajectories start from the ACDM fixed point then r increases and s decreases. 
For smaller value of interaction parameter a, the distance of {so,r } from the location of 
ACDM fixed point becomes shorter. Also the larger values of c and n result the shorter 
distance from the ACDM fixed point. The behavior of interacting polytrop ic g as in s — r 
plane is similar with the new holographic dark energy model (see Fig. (4) of [56)). For both 
models, by expanding the universe, r increases and s decreases. Finally we studied the 
w — w' analysis for this model. We showed that the evolutionary trajectories in w — w' plane 
is dependent on the parameters of the model and also the interaction parameter a. While 
the universe expands, the trajectories starts from the location of ACDM fixed point (i.e. 
{w\ = —l,w' A = 0}). Hence the polytropic gas model mimics the standard ACDM fixed 
point at the early time. The agegraphic dark energy model also mimics the ACDM model 



at the early time 40|, |4l|. At future the high-precision SNAP-type experiment can be useful 



to determine the statefinder parameters precisely and consequently single out the right dark 
energy models. 
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FIG. 1: The evolutionary trajectories for interacting polytropic gas model in s—r plane for different 
values of interaction parameter a. The black curve indicates the non-interacting case and the blue 
and red curves represent a = 0.1, 0.2, respectively. The circles on the curves show the present value 
of the statefinder pair {so, ^o}. The star symbol is related to the location of standard AC DM model 
in s — r plane. The parameters of the model are chosen as c = 2, n = 4. 
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FIG. 2: The evolutionary trajectories for polytropic gas model in s — r plane for different illustrative 
values of parameters c and n. Here we choose the interaction parameter as a = 0. In left panel 
the parameter n is fixed and the parameter c is varied as c = 2( black curve ), c = 3 ( blue curve ), 
c = 4 ( red curve ). In right panel the parameter c is fixed and the parameter n is varied as n = 2 
( black curve ), n = 4 ( blue curve ) and n = 6 ( red curve ). The circles on the curves show the 
present value of the statefinder pair {so, ^o}. The star symbol is related to the location of standard 
ACDM model in s — r plane. 
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FIG. 3: The evolutionary trajectories for interacting polytropic gas model in w\ — w'^ plane for 
different values of interaction parameter a. The black curve shows the non-interacting case and 
the blue and red curves indicate a = 0.1,0.2, respectively. The colored circles on the curves show 
the present value — iv'®. The parameters of the model are chosen as c = 2, n = 4. 
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FIG. 4: The evolutionary trajectories for polytropic gas model in u>a — w'\ plane for different 
illustrative values of parameters c and n. Here the interaction parameter is chosen as a = 0. In 
left panel the parameter n is fixed and the parameter c is varied as c = 2( black curve ), c = 3 ( 
blue curve ), c = 4 ( red curve ). In right panel the parameter c is fixed and the parameter n is 
varied as n = 2 ( black curve ), n = 4 ( blue curve ) and n = 6 ( red curve ). 
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